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' We investigate numerically, using the bond-fluctuation model, the adsorp- 

<n : 

tion of many random AB-copolymers with excluded volume interactions at 
' the interface between two solvents. We find two regimes, controlled by the 

! 

, the interface extend parallel to it, while in the second regime they extend 

c : 

Q , perpendicular to it. The density at the interface and the density in the bulk 

O ■ 



X 



depend differently on the total number of copolymers: In the first regime the 
density at the interface increases more rapidly then in the bulk, whereas the 
opposite is true in the second regime. 



PACS numbers: 61.25.Hq, 83.70.Hq 



1 



I. Introduction 

Recently copolymers at selective interfaces have received much attention [1- 
6]. Consider as an example a block copolymer consisting of a hydrophilic and a 
hydrophobic part and situated next to an oil-water interface. The difference in 
the monomers' selectivity (which we denote by kT\) favors the localization of the 
copolymer at the interface, with each block in its favorable solvent. However, for 
random copolymers frustrated situations may arise, since the chain's connectivity 
forces some monomers to stay in their unfavorable solvent. Garel et al have 
studied the localization transition of an ideal random chain at a selective interface. 
In previous publications ]7|, f| we have extended this approach to single chain under 
good solvent conditions, and we have shown that for the adsorption of random 
copolymers a simple scaling picture works very well. The scaling picture is based 
on the fact that the static properties of a single adsorbed chain consisting of N 
monomers may be understood in terms of blobs. A blob is made of g monomers; 
hence in it one of the monomer species is on the average in excess by g 1 ^ 2 . The blob 
stays in its preferred solvent as long as its total interface selectivity, ~ g^-^xksT, 
counterbalances its translational free energy, k^T . This leads to g\ 2 = C, where C 
is model dependent numerical constant. For our lattice simulation model C turns 
out to be larger than 10. Usually this constant C (which does not play any role in 
a scaling argumentation) is set to unity. It follows that the number of blobs in the 
chain, N/g, equals Ny 2 /C. In fact Nx 2 turns out to be the scaling variable of the 
problem, as we have confirmed through Monte Carlo simulations |7|, In Ref. 
we extended our analysis to asymmetric interface potentials, so that the chain as a 
whole prefers one solvent side. This problem leads to two new critical exponents (as 
predicted from scaling); we succeeded in evaluating these exponents based on our 
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simulation data 0. 

The present article is devoted to copolymer systems consisting of many chains, 
so that the polymer bulk density p& is significant. As we proceed to show, due to 
their interaction with the interface and depending on pb the copolymers display (at 
least) two distinct regimes. At very low pb the interface is only partially covered 
with adsorbed blobs; the shape of the adsorbed chains is rather flat, since their 
extension is larger parallel to the interface than perpendicular to it. In fact, the 
extension of isolated chains perpendicular to the interface does not depend on N, 
but only on the interface selectivity x [0> HI- As p& increases, the adsorbed, flat 
chains influence each other through excluded volume interactions. The chains at 
the interface form thus a two-dimensional semi-dilute solution. Increasing pb further 
leads to an interface completely saturated with blobs; thus the interface density p s 
depends in a complex manner on pb. If the adsorption is strong enough, i.e. if 
the free energy of adsorption per chain is much larger than kT, the surface can 
become saturated for values of pb for which the volume phase is still highly diluted. 
Then the adsorbed chains form large loops, resulting in an adsorption layer width 
which is of the order of the radius of gyration of the polymers in the bulk. This 



effect is well-known for homopolymer adsorption [|TJ|, n|. However, increasing pb 
beyond the saturation value a different behavior emerges: The chains extend in 
the direction perpendicular to the interface, since the loops of the adsorbed chains 
begin to stretch in a hairpin-like fashion. Such a brush-like regime for multiblock- 



copolymers at selective interfaces was recently predicted by Leclerc and Daoud |L6|. 
The reason for this behavior is that the majority blobs (or blocks in Ref.||16||) can 
be squeezed without much loss of interface energy, since only the alternation of the 
two blob types across the interface fixes the chain. In this way more polymer chains 
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can be adsorbed at the interface. 

Our simulation results confirm this picture qualitatively. Distinct from the be- 
havior of adsorbed homopolymers, we report here for pb larger than a characteristic 
value pi the stretching of the adsorbed copolymers in the direction perpendicular to 
the interface and their contraction parallel to it. 



II. Simulation algorithm 

Our Monte Carlo simulations for copolymers were performed using the bond- 



fluctuation method (BFM) JT0|, 11]. The BFM is a lattice algorithm where each 
monomer is represented by a lattice cell. Thus on a three-dimensional simple cubic 
lattice each monomer occupies eight neighboring lattice sites. The length of a bond 
connecting two neighboring monomers fluctuates between 2 and \/lQ lattice spac- 



ings |TIJ. Self-avoidance (the excluded volume interaction) is satisfied by not 
allowing any two monomers to have a lattice site in common. To avoid bond-crossing 



the allowed bonds are restricted to a set of 108 vectors 0, [LI], [T^l . Here we study 
the behavior of n random copolymers, each of length N, all placed in a box of size 
L x L x H, with periodic boundary conditions in the x- and y-directions and two 
impenetrable surfaces at z = and z = H. The copolymers are random; for each 
of them separately the iV- monomer sequence consists of randomly chosen A- and 
5-monomers. For adsorption we assume a symmetrical situation: the interaction 
energy of each monomer with its unfavorable solvent is x^bT and with its favorable 
solvent is zero. We let the solvent below the interface (z < H/2) favor A-type, 
and the solvent above the interface (z > H/2 + 1) B-type monomers. Note that 
the interface is thus at zq = (H + l)/2. In the Monte Carlo algorithm the chains 
move by position changes of their monomers, which attempt nearest neighbor steps 
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on the underlying cubic lattice. A move is taken into consideration only if it satis- 
fies the requirements of self-avoidance and of non-breaking of bonds. Furthermore, 
energetically unfavorable moves are statistically permitted according to the usual 
prescription involving the Boltzmann factor. 

III. Simulation Results 

We study the density dependence of the adsorption properties by changing n, 
the number of polymers in the LxLxH box. Here we take L = 50 and H = 100 and 
focus on the results obtained using copolymers of length N = 64 with a monomer- 
solvent interaction parameter of x — 3.15. Results for other copolymer lengths and 
for other \ parameters will be mentioned when appropriate. Note that in the single 
chain case (n = 1) the parameters N = 64 and x = 3.15 let the system be located 
in the well adsorbed scaling regime, see Ref. [0, ||. 

An initial configuration is generated by randomly placing the first monomer of 
each polymer in the system and then randomly adding the subsequent monomers, 
such that self-avoidance and the restrictions on the bonds are obeyed. The energetic 
aspects of the interaction with the solvents are then taken care of by the usual 
Boltzmann factor; the monomer-monomer interaction is only accounted for through 
the excluded volume aspect. This means that both solvents are good for both 
species. We established numerically that the relaxation time (determined using 
the autocorrelation function of the radius of gyration R g and of its z-component 
Rg± 0) is around 50,000 Monte Carlo steps (MCS), where a MCS consists of nN 
move attempts; we thus view the copolymers as having reached equilibrium after 
200,000 MCS. Averages are then calculated from the configurations obtained in the 
subsequent 200,000 MCS. In order to improve the statistics we average over results 
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from 8 independent runs. 

In Figure 1(a) we plot the densities of the A-type monomers (pa) and of the 
B-type monomers (pb) as a function of the height z. We pause to make clear how 
these densities are normalized. In the lattice, bond-fluctuation model used here, the 
maximal number of monomers which can be accommodated in the given volume is 
^L 2 H, since due to the excluded volume restrictions each monomer blocks 8 lattice 
sites. In our system containing n chains of N monomers each the total number of 
monomers is nN; hence the average density is p = 8Nn/(L 2 H). Now the densities 
Pa and ps satisfy the relation 

-EtT,M z ) + Pb(*)) = P ■ (!) 

" z 

The number of polymers in Fig. 1(a) is n = 50. Comparing with the results for 
single chain, Fig.l of Ref.lTJ, we find that now the monomer densities are still nonzero 
quite far away from the interface (i.e. from z = 10 to z = 30 and from z = 70 to 
z = 90 in Fig. 1(a)). Furthermore, at such distances pa and ps are equal and 
are independent of z. In this range we identify this constant with the bulk density. 
When approaching the impenetrable boundaries at z = and at z = 100, pa and pB 
drop to zero. In the following we will not consider the range from z = to z — 10 
and from z = 90 to z = 100 any further. Close to the interface pa and ps peak 
on their favorable side. The densities decay smoothly on their favorable side and 
sharply across the selective interface, so that their values on the unfavorable side 
near the interface lie below the bulk density. This differs from our findings for single 
chain, for which we found that the bulk density is zero and that densities close to 
the interface on the unfavorable side display a secondary peak O. In Fig. 1(b) we 
show pa + Pb as a function of z. The result is a symmetric peak centered at the 
interface and superimposed on a bulk density background. 
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In Fig. 2(a) we plot the bulk density pb against the average density p, as given 
by Eq.flU). We see that pb becomes extremely low for large x an d small p (but does 
not disappear completely). The density pb increases linearly with p for larger p. 
Fig. 2(b) shows the density at the interface p s = (pa(H/2) + pn{H/2) + pa(H/2 + 
1) + Pb(H/2 + l))/2 as a function of p. For small p, p s increases rapidly, while 
for larger p it increases more slowly, the cross-over region being around p ~ 0.02. 
From the Figure the almost linear increase of p s in the range of large p is also clear. 
Fig. 2(c) shows a plot of p s versus pb, where again the two regimes can be seen. To 
understand these regimes we plot in Fig. 2(d) p s — p b against p. In terms of Fig. 1(b) 
p s — pb is in fact the peak's height relative to the background. We see from Fig. 2(d) 
that p s — pb attains its maximum around p c ~ 0.02. For p < p c , p s — p b increases 
with increasing p while for p > p c it decreases with increasing p. This means that 
in the low p regime (p < p c ) adding more polymers to the system leads mainly to 
an increase in density at the interface, whereas for higher p values (p > p c ) adding 
more polymers to the system leads to an overall increase in density in the bulk. In 
the second regime both p s and pb increase almost linearly with p, whereas p s — pb 
decreases roughly linearly with increasing p. Extrapolating the linear dependence 
in the second regime of Fig. 2(d) to p s — p b = leads to p ~ 0.95. We note that a 
third regime may exist at still higher p values, when the copolymer concentration in 
the bulk reaches the semi-dilute range, but investigations in exploring rather dense 
systems are beyond the scope of the present paper. 

For the record, we like to point out that varying the copolymers' length, using 
N = Q4 and iV = 128, we obtained plots almost identical to those shown in Fig. 2. 
For different x parameters (x = 2.20, 3.15 and 4.20), we find that pb is independent 
of x while p s increases with increasing x- 
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The difference between pa and pB provides a means to quantify the interfacial 
selectivity. Fig. 3(a) displays the data for n = 50. We plot Pa-b = Pa — Pb f° r 
z < H/2 and pa-b = Pb — Pa for z > H/2 + 1; this leads to a symmetric peak. We 
find in this density regime (p ~ 0.1) that the width of the peak is independent of N 
but is controlled by x- Fig. 3(b) shows how the width of the peak in Fig. 3(a) varies 
as a function of p for different x- Here the width is taken as full width at half-height 
(FWHH), i.e. the difference Z2 — zi, with pa-b(zi) = Pa-b{z-i) = \pa-b{zq)i where 
Pa-b(zq) is the density at the interface (the peak's height in Fig. 3(a)). Since our 
model is based on a lattice, we had to interpolate between the discrete values, in 
order to obtain the (non-integer) FWHH. We have also evaluated numerically the 
second moment of the distribution of Fig. 3(a) and found indications that the second 
moment may diverge. In Fig. 3(c) we show the peak's height pa-b as a function of 
p and we find that the relationship is linear above p c . 

Keeping in mind that the interfacial selectivity affects mostly the chains close to 
the interface, i.e. the adsorption layer, see Fig. 3(b), we turn now to a comparative 
study of the polymers' behavior around the interface and also away from it. For 
this we compute R g ^ and R g \\, the z- and the xy-components of the radius of 
gyration of each copolymer. The center of mass of the fcth polymer is Rqm = 
J2iLi Ri®/N, where = (x\ k \ y\ k \ z\ k ^) are the coordinates of the ith monomer 
within the fcth chain. For the fcth chain, R g k [ = (E^Ii^P ~ z^l) 2 /N) 1 / 2 , and 
4? = (I ZliK^ - xci) 2 + {y? ~ y { c k) M ) 2 ]/N) 1/2 . Now we have to specify which 
copolymers belong to the space around the interface. For this we use a a-criterion: 
If the z-component of the polymer's center of mass is within a distance a from the 
interface, i.e. if \z^ M — z \ < c, we view the polymer as being near the interface, 
otherwise as being far from it (here again we disregard the far-off regions close to the 
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fixed boundaries). Now, evidently 2a should be taken larger than the peak's FWHH. 
Since the FWHH, Fig. 3(b), depends on p, we choose a = 10 in what follows (we 
checked that the choice of a does not change the features reported below, provided 

(k) (k) 

that 2o" is reasonably larger than the FWHH). We now average R>gl and R$ for the 
polymers near the interface to obtain R g ± and R g \\\ we do the same for the polymers 
far from the interface, which leads to R g ± and R g \\. In Fig. 4(a) and Fig. 4(b) we 
present plots of R g ± and R g \\ and of R g ± and R g \\ as functions of p. 

Let us first consider the bulk phase, away from the interface. As can be readily 
inferred from Fig. 4(b), R g ± is close to R g » for all p, which means that the bulk 
copolymers' shape is not affected by the interface and that it is isotropic in space. 
Note that the large fluctuations of the data points for small p values arise from the 
fact that almost no chains are located in the bulk until the interface is saturated, 
see also Fig. 2(a). 

The situation is completely different for polymers near the interface, Fig. 4(a). 
For very small values of p, R g ± is considerably smaller than R g \\; hence adsorbed 
chains tend to be rather flat. Fig. 4(a) shows that R g » decreases monotonically with 
increasing p, whereas the opposite is true for R g ±. This is analogous to the behavior 
of adsorbed homopolymers, as discussed by Bouchaud and Daoud |T4]| . Comparing 
Fig. 4(a) with Fig. 2(c) one notes that the chain's extension is, in contrast to the 
surface density, a smooth function of p also at the saturation value p c . However, 
increasing the density beyond the value p* ~ 0.05, where R g u = R g ±, the extension 
perpendicular to the interface gets to be larger than the parallel extension. This 



is in accordance with Ref. flB |, which predicts (as discussed in the Introduction) 
that above a characteristic bulk concentration the loops of the adsorbed chains will 
stretch in the direction perpendicular to the interface, in a hairpin-like fashion. 
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Note that also the chain's extension parallel to the interface, R g », decreases below 
its bulk value R g \\ (compare Figs. 4(a) and (b)), i.e. the adsorbed chains begin to get 
squeezed parallel to the interface. Thus varying the bulk concentration only by a 
small amount (0 < pb < 0.1) the chain's geometry changes from a flat, pancake-like 
shape into a brush-like assembly of stretched loops. This picture is also supported 
by the behavior of p s clS db function of p in the saturated surface regime (p > p c ). As 
can be inferred from Fig. 2b, the surface concentration increases even beyond the sat- 
uration threshold p c (which is given in Figs. 2 already by the second numerical point). 

IV. Conclusions 

We have investigated the behavior of many random copolymers in the presence 
of a selective interface. Previous studies || ||, |7|. || showed that a simple scaling 
picture works very well for random copolymer adsorption. In this paper we have 
extended our previous work to consider the effects of sizeable chain concentrations. 
We found that there are at least two regimes controlled by the polymer density, as 
can be inferred from Figs. 2 and Fig. 4a. 

In the very low density regime (except for the adsorption mechanism) copoly- 
mers at interfaces behave similarly to homopolymers at surfaces. Starting from a 
single chain, which in the adsorbed state is flat, an increase in the chains' density 
leads to their crowding at the interface. As usual for polymer adsorption, the inter- 
face may be fully covered by chains even when the bulk density is still highly diluted 
(compare the surface density p s and the bulk density pb for the second numerical 
point in Fig. 2(c)). This can be easily understood from the fact that the adsorption 
energy per chain is in most huge quantity compared to the translational free 

energy per chain, i.e. to kT. 
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In the second regime, however, where the interface is already covered, the only 
way of adding more chains to it is to squeeze the chains in a brush-like fashion. As 
a consequence, R g ±, the radius of gyration of the adsorbed chains perpendicular to 
the interface exceeds the average value R g ± in the bulk. A means of picturing this 
situation is a brush-like assembly of stretched loops [16| . 
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Figure Captions 

Figure 1(a): Densities of the A-monomers (pa) (diamonds) and of the B-monomers 
(Pb) (crosses) at the height z. Here the number of polymers is n = 50, the length of 
all polymers is N = 64, and the monomer-solvent interaction parameter is x — 3.15. 

Figure 1(b): Display of \pa + Pb\ vs. z, see Fig. 1(a) for definitions. 

Figure 2(a): The bulk density (p b ) vs. the average density (p), see text for details. 

Figure 2(b): The density at the interface (p s ) vs. p. 

Figure 2(c): p s vs. p b . 

Figure 2(d): p s - p b vs. p. 

Figure 3(a): pa-b (the absolute value of pa — Pb) vs. the height z. Here n = 50. 

Figure 3(b): The FWHH of Fig.3(a) vs. p for three different x parameters, x = 2-20 
(diamonds), x = 3.15 (crosses), x — 4.20 (squares). 

Figure 3(c): Pa-b(zq) vs. p, with z being the interface's location. 

Figure 4(a): Radius of gyration for polymers near the interface: R g j_ (diamonds) 
and R g \\ (crosses) vs. p. 
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Figure 4(b): Radius of gyration for polymers away from the interface: R g ± 
monds) and i? s || (crosses) vs. p. 
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